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ABSTRACT 

Context. A subclass of the upper main sequence chemically peculiar stars, mercury-manganese (HgMn) stars were traditionally 
considered to be non-magnetic, showing no evidence for variability in spectral line profiles. However, recent discoveries of chemical 
inhomogeneities on their surfaces call for revision of this question. In particular, spectroscopic time-series of AR Aur, a And, and five 
other HgMn stars indicate the presence of chemical spots. At the same time, no signatures of global magnetic fields were detected. 
Aims. In order to understand the physical mechanism that causes the formation of chemical spots in HgMn stars and to gain insight 
into the potential magnetic field properties at their surfaces, we performed a highly sensitive search for magnetic fields on a large set 
of HgMn stars. 

Methods. With the aid of a new polarimeter attached to the HARPS spectrometer at the ESO 3. 6m- telescope, we obtained high-quality 
circular polarization spectra of 41 single and double HgMn stars. Using a multi-line analysis technique on each star, we co-added 
information from hundreds of spectral lines resulting in significantly greater sensitivity to the presence of magnetic fields, including 
very weak fields. 

Results. For the 47 individual objects studied, including 6 components of SB2 systems, we do not detect any magnetic fields at greater 
than the 3<r level. The lack of detection in the circular polarization profiles indicates that if strong fields are present on these stars, 
they must have complex surface topologies. For simple global fields, our detection limits imply upper limits to the fields present of 
2-10 Gauss in the best cases. 

Conclusions. We conclude that HgMn stars lack large-scale magnetic fields, typical for spotted magnetic Ap stars, sufficient to form 
and sustain the chemical spots observed on HgMn stars. Our study confirms that in addition to magnetically altered atomic diffusion, 
there exists another differentiation mechanism operating in the atmospheres of late-B main sequence stars which can compositional 
inhomogeneities on their surfaces. 
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1. Introduction 

Mercury-manganese (HgMn) stars form a subclass of the upper 
main sequence chemically peculiar (CP) stars, showing a no- 
table overabundance of Hg, Mn, Y, Sr and other, mostly heavy, 
chemical elements with respect to the solar chemical compo- 
sition. HgMn stars are frequently found in binaries and lie on 
the H-R diagram between the early-A and l ate-B spectral type s, 
which corresponds to r eff = 9500-16000 K (iDworetskvi 1 1 9931) . 

While most other, normal stars in this spectral range are rapi d 
rotators with v e sin i = 200 - 300 km s -1 dAbt & MorrellL fl995). 
HgMn stars are typically slow rotators, making them ideal late- 
B targets for detailed abundance analyses. The HgMn stars were 
for a long time considered to have chemically homogeneous at- 
mospheres, with no obvious vertical chemical abundance strati- 
fication or horizontal concentrations of some elements into sur- 
face spots of high elemental abundan ce such as those in mag- 
netic Ap stars of similar temperatures (Kochukhov, 2004). This 
view was historically supported by the lack of a definite detec- 
tion o f line profile variabilit y in any HgMn star. However, re- 
cently Adelman et al] (120021) reported variability in the spectral 

* Based on observations collected at the European Southern 
Observatory, Chile (ESO programs 083.D-1000, 084.D-0338, 085.D- 
0296). 



line of Hg n A 3984 A based on high-resolution and high signal- 
to-noise (S/N) ratio spectroscopic data obtained for the brightest 
HgMn star: a And. They attributed the variability in the Hg n line 
to horizontal inhomogeneities of the mercury abundance across 
the stellar surface and reconstructed a surface map of Hg with 
the help of Doppler imaging. 

Later, IKochukh ov et al. I J2001 found two other HgMn stars 
(HR 1185 and HR 8723) possessing spotted structure in their 
abundance of Hg. The fourth star, the eclipsing HgMn binary 
AR Aur studied bv lHubrig et al.l d2006l) . was reported to exhibit 
line profile variability of Hg and also such chemical elements as 
Y, Zr, Pt, an d Sr. These results w ere confirmed by the indepen- 
dent study of iFolsom et alJ (l2010h . 

An addition a l une xpected discovery was made by 
IKochukhov et all (120071) for the HgMn star a And. From 
the analysis of the Hg 11 3984 A line observed over a timespan of 
7 years they inferred surface maps which display morphological 
changes of spots on a time scale of only a few years. In 
comparison to well-studied magnetic CP stars, this behavior of 
the surface abundance distribution is quite surprising because 
all other early-type spotted stars show spot configur ations that 
are stable over tens of years (e.g jAdelman et all 12001 ). 
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Finally, Briqu et et al.l (1201 Oh investigated time-series spectra 
of three more HgMn stars and concluded that all of them dis- 
play variability in their line profiles. The star HD 1 1753 showed 
variability in Y, Sr and Ti; while HD 53244 and HD 221507 
both show variable profiles of their Mn and Hg lines. In addi- 
tion, HD 221507 displays variations in its Y lines. In summary, 
these recent studies show there are at least 7 spotted HgMn stars 
and possibly many more. It is still not clear how typical this be- 
havior is for this class of CP stars. Furthermore, it is not clear 
what physical processes are responsible for the formation of 
these structures on the stellar surface. 

It is generally believed that a strong magnetic field is a nec- 
essary ingredient for creating inhomogeneities in the stellar at- 
mosphere, leading, for e xample, to the formation of temperature 
or ch emical spots (e.g., Jeffers & Donati, 2008; Liiftingeret al., 
l20loh . Within this framework, a number of spectropolarimet- 
ric observations have been carried out on HgMn stars in order 
to confirm or disprove the presence of magnetic fields on these 
stars which could be responsible for the formation of chemical 
spots. The first systematic attempt to detect magnetic signatures 
in high-resol ution circularly pol arized spectra of HgMn stars 
was made bv lShorlin et al.l y002). This study included a sample 
of 10 stars which were analysed at a resolution of R = 35 000, 
yielding upper limits for the magnetic field strength of 29 to 
100 G. To date, this is the only extensive magnetic survey of 
HgMn stars. Several other studies searched fo r magnetic f i elds in 
individual spotted HgMn stars. For instance. IWade et all {2006) 
searched for magnetic a field in a And. They analysed magnetic 
field measurements obtained with three different polarimeters 
and concluded that the sta r has no longitudina l magnetic field 
stronger than about 6-19 G. lFolsom et al.l (|2010) performed sim- 
ilar analysis for AR Au r, once again findin g no field stronger 
than 20-40 G. Finally. lAuriere etaTI (Hp 10} presented an anal- 
ysis of a small sample of three bright, sharp-lined HgMn stars, 
reporting no magnetic field detections at the level of 1-3 G. 

Previous magnetic field studies of HgMn stars suffered from 
several fundamental limi tations. Firs t , with the exception of 
the fo ur stars studied by IWade et al] d2006l) and lAuriere et alj 
(1201 Oh . previous work was not particularly precise, possibly 
missing weak magnetic fields which are still capable of causing 
chemical spot formation. Second, they included a small number 
of HgMn stars, often with a strong emphasis on the most slowly 
rotating ones. Thus, the class as a whole has not been well sur- 
veyed. 

To overcome drawbacks associated with previous magnetic 
field studies on these stars, and to provide new insights into the 
spot formation physics of late-B stars, we have initiated a spec- 
tropolarimetric survey of a large sample of HgMn stars. We in- 
vestigate HgMn stars with a broad range of atmospheric param- 
eters and rotational velocities. Most of these objects have not 
been studied before with the high-resolution spectropolarimetry. 
Taking advantage of a new polarimeter, HARPSpol, attached to 
the HARPS instrument at the ESO 3.6-m telescope in La Silla, 
Chile, we were able to push the limits of magnetic field detection 
in early-type stars down to remarkably low levels. 

In Sect. 12.11 we describe our target selection proce- 
dure. Sect. 12.21 discusses our spectropolarimetric observations. 
Sect. 12.31 outlines the principles of the reduction of the spec- 
tropolarimetric data from HARPSpol. We discuss our multi-line 
magnetic field detection technique in Sect. 13.11 and present the 
results of our magnetic field measurements in Sect. 13.21 Sect. [4] 
summarizes our results and discusses them in the context of re- 
cent studies of HgMn stars. 



Table 1. The list of HgMn stars included in our survey. 



HD number 


HR number 


my 


T eB (K) 


Binarity 


HD 1909 


HR 89 


6.56 


12406 


SB2 


HD 11753 


HR558 


5.11 


10476 


SB1 


HD 27376 


HR 1347 


3.55 


12849 


SB2 


HD 28217 


HR 1402 


5.87 


13906 


SB1 


HD 29589 


HR 1484 


5.45 


14763 




HD 32964 


HR 1657 


5.10 


11119 


SB2 


HD 33647 


HR 1690 


6.67 


12440 


SB2 


HD 33904 


HR 1702 


3.28 


12759 




HD 34880 


HR 1759 


6.41 


13269 




HD 35548 


HR 1800 


6.56 


11164 


SB2 


HD 36881 


HR 1883 


5.63 


10903 


SB1 


HD 42657 


HR 2202 


6.20 


12842 




HD 53244 


HR 2657 


4.10 


13596 




HD 53929 


HR 2676 


6.09 


13908 




HD 63975 


HR 3059 


5.13 


13458 




HD 65950 




6.87 


12710 




HD 68099 


HR 3201 


6.08 


12997 




HD 70235 


HR 3273 


6.43 


12329 




HD 71066 


HR 3302 


5.62 


12131 




HD 71833 


HR 3345 


6.67 


12985 




HD 72208 


HR 3361 


6.83 


11141 


SB2 


HD 75333 


HR 3500 


5.31 


12248 




HD 78316 


HR 3623 


5.24 


13639 


SB2 


HD 90264 


HR 4089 


4.95 


15121 


SB2 


HD 101189 


HR 4487 


5.14 


11148 




HD 106625 


HR 4662 


2.59 


12002 


SB1 


HD 110073 


HR4817 


4.63 


12876 




HD 141556 


HR 5883 


3.96 


10684 


SB2 


HD 149121 


HR6158 


5.62 


11021 


SB1 


HD 158704 


HR 6520 


6.06 


13378 


SB2 


HD 165493 


HR 6759 


6.15 


14375 


SB2 


HD 175640 


HR7143 


6.20 


12075 




HD 178065 


HR 7245 


6.56 


12348 


SB1 


HD 179761 


HR 7287 


5.14 


13010 




HD 186122 


HR 7493 


6.33 


12901 




HD 191110 


HR 7694 


6.18 


12107 


SB2 


HD 193452 


HR7775 


6.10 


10881 


SB1 


HD 194783 


HR7817 


6.08 


13803 




HD 202671 


HR8137 


5.39 


13696 




HD 211838 


HR8512 


5.35 


12593 


SB1 


HD 221507 


HR 8937 


4.37 


12476 





2. Observations 

2.1. Target selection 

To compile our target list we used the catalogue of A p, HgMn, 
and Am stars published bv lRenson & Manfroidl d2009l) . In order 
to optimize our observations for the detection of weak magnetic 
fields and to cover a statistically large sample of HgMn stars, 
we put several constraints on the stellar parameters. All objects 
with v e sin / ^70 kms~' were excluded because magnetic field 
measurements are considerably less accurate for rapid rotators 
in comparison to the moderately and, especially, slowly rotating 
stars. We also preferentially selected brighter stars. Aiming to 
reach a signal-to-noise ratio (S /N) of « 300, we selected only 
targets brighter than my — 7. A total of 45 HgMn stars were 
found that satisfied these criteria. Our three observing runs cov- 
ered 41 of them. Among these stars 6 are spectroscopic binary 
systems, allowing magnetic field measurements on both compo- 
nents. 

In the Table Q] we present the list of observed stars, giv- 
ing their HD and HR numbers, visual magnitude, effective tem- 
perature and binary status. The effective temperature of all ob- 
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jects was estimated using the iMoon & Dworetskvl (1 19851) cal- 
ibration of the Stromgren photo metric paramete r s as i mple- 
mented in the TempLogG code (Kup ka & Brunttl E001). The 
Stromgren photometry was obtained from the SimbacQ database. 
Information ab out binarity was extracted from the catalogue of 
Pour baix et alj d2004l) and complemented by the results of sev- 
eral recent studies of individual binary systems. For a few stars 
we used our LSD line profiles as described below to detect bina- 
rity. 

2.2. Spectropolarimetric observations 

Our observatio ns were obtained u sing the newly built polarime- 
ter HAR PSpol dSnik et all l2008h attached to the HARPS spec- 
trometer (iMavor et al.L 120031) at the ESO 3.6-m telescope. With 
a resolving power of R = 115 000, this instrument is the high- 
est resolution spectropolarimeter available to the astronomical 
community. The compact optical design of the polarimeter al- 
lows it to be mounted before the fiber entrance at the Cassegrain 
focus of the telescope. This positioning of the polarimeter min- 
imizes the instrumental polarization which is usually caused by 
the oblique reflections in the light path. HARPSpol consists of 
two independent polarimeters for circular and linear polariza- 
tion measurements. Thus, it is possible to make observations in 
all four Stokes parameters. Each polarimeter consists of a super- 
achromatic retarder plate (a half-wave or a quarter-wave plate) 
and a beam splitter based on a Foster prism. This design makes 
the polarimeter fully achromatic. The two polarized light spectra 
are recorded simultaneously through the two HARPS entrance 
fibers. 

Our observations were carried out between May-June of 
2009, January of 2010, and April-May of 2010. The date each 
of our target stars was observed is recorded in Table 2. In to- 
tal, we have observed HgMn targets on 26 nights, which were 
shared with other spectropolarimetric programs. We obtained 
high-resolution spectra with a typical S /N of 150^-00. The 
spectra were recorded by a mosaic of two 2K x 4K CCDs, pro- 
viding 45 polarimetric echelle orders on the blue CCD and 26 
on the red one. For each night we have acquired a standard set of 
calibration images: 20 bias exposures, 20 flat fields and 2 ThAr 
frames. Flat field and ThAr images were exposed with the po- 
larimeter in the circular polarization mode. 

All HgMn stars were observed in circular polarization, cov- 
ering the wavelength range of 3780-6913 A with a small gap at 
5259-5337 A. Each observation of an individual star was split 
into four sub-exposures corresponding to differing positions of 
the quarter-wave plate: 45°, 135°, 225° and 315° relative to the 
optical axis of the beam-splitter. The length of individual sub- 
exposures ranged between 160 and 350 s. In a few cases stars 
were observed with only 2 sub-exposures. The majority of our 
targets were observed only once; however, 10 stars, mainly spec- 
troscopic binaries, were observed several times. For HD 11753 
and HD 32964 we obtained complete coverage of the rotational 
and orbital periods respectively. Analysis of these HgMn stars 
will be presented in future publications. 

2.3. Reduction of the data 

The reduction of our HARPS pol observations was perf ormed 
with the REDUCE package of jPiskunov & Valentil d2002l) . This 
set of IDL routines applies a standard sequence of reduction 
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and calibration procedures to the cross-dispersed echelle spec- 
tra. Bias images are averaged and subtracted from the average 
flat and science images. Spectral orders are located using the av- 
erage flat field with the help of a cluster analysis method. The flat 
field is normalized and is used to correct the pixel-to-pixel sensi- 
tivity variations in the science images. After the removal of scat- 
tered light, the science spectra ar e extracted using the op timal 
extraction algorithm described bv lPiskunov & Valentil (120021) . 

The intrinsic long-ter m stability of the H ARPS spectrometer 
is of the order of 1 m s -1 (M ayor et al.ll2003l) . We do not expect 
the polarimeter to adversely affect this stability because it does 
not significantl y modify the see ing-limited illumination of the 
HARPS fibers (Sni k et aTll2008h . Since we do not require an ex- 
tremely high velocity accuracy for our observations, it was suffi- 
cient to use only one nightly ThAr spectrum for the wavelength 
calibration. Using « 700-900 ThAr lines in a 2-D wavelength 
calibration routine, we obtained an internal wavelength calibra- 
tion accuracy of 18-21 m s . In the final step before the cal- 
culation of the Stokes parameter spectra we performed contin- 
uum normalization, dividing each spectrum by a smooth, slowly 
varying function. This function was obtained by fitting the upper 
envelope of the blaze-corrected merged spectrum. 

For the calculation o f the circular polarizat ion we use the ra- 
tio method described by Bag nulo et al.l (120091) . The method re- 
duces the spurious polarization effects by an appropriate com- 
bination of the two physical beams recorded for four different 
retarder plate positions. Along with the circular polarization we 
also derive a diagnostic null spectrum. It is obtained by com- 
bining individual sub-exposures destructively, thus canceling the 
stellar polarization signal and showing the residual instrumen- 
tal polarization remaining after application of the ratio method. 
Below we use the null spectrum in the same analysis steps as the 
Stokes V spectrum, thus providing a realistic estimate of possi- 
ble errors. 

Using the eight spectra in the four individual sub-exposures 
of each star provides a convenient means to detect and remove 
cosmic ray hits that otherwise seriously distort the final Stokes 
spectra. Affected pixels are identified by their large deviation 
from the median and are substituted by the latter. Since our ex- 
posure times were relatively short, only 1-2 pixels required cor- 
rection in each echelle order. 



3. Data analysis and results 

3.1. LSD analysis 

A major difficulty when searching for weak (below 100 G) mag- 
netic fields in stars is the weakness of expected polarimetric 
signal. As one can see from Fig. [TJ it is difficult to conclude 
anything about the presence or absence of the magnetic sig- 
nal in individual spectral lines even with our best Stokes V 
spectra. To alleviate this problem, a line addition technique, 
ca lled Least-Squares D econvolution (hereafter LSD) developed 
bv lDonati et al.l (1 1 997b . is commonly used to detect weak stellar 
magnetic fields. This technique has proven to be a very effective 
tool, allowing extraction of high-precision Stokes / and V pro- 
files by combining the information from all available metal lines 
in the spectrum. The main assumption of LSD is that all spectral 
lines are identical in shape and can be represented by a scaled av- 
erage profile. Using this model one can reconstruct the average 
profile from an observed intensity or polarization spectrum, ob- 
taining an increase of the signal-to-noise ratio by a factor of up to 
30—40. Recent publications have demonstrated that with the help 
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Fig. 1. The HARPSpol spectra of HD 71066 in the wavelength region of 5220-5240 A. Spectra are displayed from bottom to top 
as follows: Stokes /, Stokes V and the null spectrum. Note, that for displaying purposes Stokes V and the null spectrum are shifted 
upwards by 1.12 and 1 .25, respectively, and are multiplied by 5. 



of the LSD techniq ue, it is possible to measure magnetic fie lds 
weaker than 1 G (lAuriere et a l.. 2009; Lignieres et all 120091) . 

The LSD analysis requires prior knowledge of the line posi- 
tions and st rengths. We have com piled a number of line lists us- 
ing VALD dPiskunov et al.lll995l) together with a 1-D, LTE stel- 
lar atmosphere models with effective temperatures ranging from 
10000 K to 18000 K in 500 K steps. T he atmosphere model s 
were calculated with the LLmodels code (IShulvak et all [2004), 
assuming solar abundances. When doing the spectrum synthe- 
sis to estimate central line depths in order to select lines to in- 
clude in our line lists, the abundances of some chemical elements 
were changed with respect to the Sun to mimic the abundance 
pattern typical of an HgMn star, e.g. [Fe] = +0.2, [Cr] = +0.5, 
[Mg] = -0.2, [Ti] = +0.5. 

The effective temperature of each star, derived from the 
Stromgren photometric colors, was rounded to the nearest value 
in the model grid. A separate line mask was then created for each 
stare. In order to exclude very weak lines from the calculation of 
LSD profiles, we required that the central depth of each line be 
> 10% of the local continuum in order for the line to be included 
in the list. For the coolest star in our sample (T e g « 10500 K) 
this gave us 753 lines for the calculation of the LSD profiles, 
while for the hottest star (T ef( « 14000 K) we used 344 lines. 
The resulting gain factor ranges from 5 to 17, depending on the 
S /N of observations and the T e ff of the star. 

For Stokes /, the line weights employed in the LSD analysis 
are assumed to be the central line depths provided by VALD. For 
the Stokes V spectrum, we computed line weights using 

w = dz^-, (1) 

where d is the central line depth, z is the effective Lande fac- 
tor, A is the laboratory wavelength of the respective line and Aq 
is a normalization wavelength, which we assumed to be equal 
to 4800 A. In a few cases when Lande factors of the upper 
and/or lower atomic levels were missing in VALD, we used 
the LS-coupling scheme to compute them. In all other cases 
Lande factors come from the quantum-mechanical calculations 
of Kuruc20. This way of assigning weights is considered stan- 
dard (|DonSiitall[l993). 

2 http://kurucz.harvard.edu/linelists.html 



We derived the Stokes /, V and null LSD profiles using the 
code written by one of us (Kochukhov et al., submitted). All 
profiles were reconstructed using a velocity grid with a bin- 
size of 0.8 kms _I , corresponding to the average pixel scale of 
HARPSpol spectra. The formal uncertainties of the LSD profiles 
were obtained by propagating the variance provided by the spec- 
tral reduction code for individual pixels. To account for possible 
under- or over-estimation of the error bars and for the eventual 
failure of the LSD assumptions, we scaled the uncertainties of 
the calculated LSD profiles to achieve redu ced chi-square val ues 
of 1 .0, in the same manner as described bv lWade et al.l d2000t) . In 
Fig. [2] we show an example of the LSD spectra obtained for the 
sharp-lined HgMn star HD 71066. This is the object for which 
we have achieved our lowest uncertainty for the magnetic field 
diagnostics. 




-20 20 40 60 

Av, km s" 1 



Fig. 2. The resulting LSD profiles for HD 71066. For displaying 
purposes LSD V and null profiles are shifted upwards by 1.08 
and 1.17, respectively, and are multiplied by 25. 
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3.2. Magnetic field measurements 

The mean longitudinal magnetic field, (B z ), was inferred from 
LSD profiles. Prior to the magnetic field measurements, we 
renormalize each LSD profile by a constant factor in order to 
obtain the correct continuum normalization. The displacement 
of the continuum is caused by the contribution of numerous un- 
accounted for weak blends in the LSD / profile, resulting in a 
slightly lower continuum level. The scaling of LSD /, V and null 
profiles was done consistently for each star. Then renormalized 
LSD profiles were used to estimate the longitudinal magnetic 
using 

f V(v - v )dv 

(B z ) = -7.14 x lO 6 ^ , (2) 

J (1 - I)dv 

where V corresponds to the LSD circular polarization spectrum 
and / corresponds to the LSD intensity profile, v is the veloc- 
ity, and vo is the velocity of the center-of-gravity of the LSD 
Stokes / profile. The constant factor in the right hand side of 
Eq. (O includes Ao-the normalization wavelength of the LSD 
profile, equal to 4800 A, and zo - Lande factor of the LSD pro- 
file, equal to 1 according to our definition in Eq. (Q]i. The re- 
sulting values of (B z ) are in Gauss. The uncertainties in the (B z ) 
measurements are calculated by propagation of the LSD profile 
uncertainties. 

The integration limits in Eq. (|2]) were chosen on the basis of 
the LSD / profile for each star. We note that the choice of inte- 
gration limits can be quite important for the resulting (B z ) and 
its errors bar, especially for stars with v e sin i SS 25 kms -1 . The 
integration limits must be selected symmetrically with respect to 
the core of the LSD / profile. They should not exclude the wings 
of the profile but, at the same time, should not extend too far into 
the continuum which contains no useful polarization signal. We 
find that for very slowly rotating stars the sensitivity of our (B z ) 
measurements can be degraded by 5-7 G if the integration limits 
are changed by ~3 kms -1 . In this case we carefully chose the 
integration limits in order to reach the best precision of the lon- 
gitudinal field determination. In the case of faster rotators, the 
definition of the measurement window is less important, so the 
integration limits can be set with the precision of ~5 km s . 

The same measurement procedure was also applied to the 
LSD profiles calculated from the null spectrum in order to test 
whether a real polarization signature is contained in the stellar 
spectrum or if it is produced by spurious polarization. In general, 
the uncertainties of the longitudinal magnetic fields determined 
from the Stokes V and null spectra are in excellent agreement 
and in none of the cases do we derive significant longitudinal 
field from the null profile. This result confirms the quality of our 
observations and analysis and suggests the absence of noticeable 
spurious polarization. 

We have measured magnetic fields for 41 objects. Since the 
target list includes SB2 systems with the contribution of both 
stars clearly visible in the LSD profiles, we also measured mag- 
netic fields in 6 secondary components. Most of them are Am 
stars dRvabchikoval 1 1 9981) . except HD 27376 where the sec- 
ondary has Hg peculiarity. There are 68 (B z ) measurements in 
total, as some of the stars were observed more than once. 

For all stars, our measured (B z ) values are below the 3<x 
level. Similarly we do not detect (B z ) using the null spectra. The 
best precision we achieved for our longitudinal magnetic field 
measurements was for HD 71066, resulting in an uncertainty of 
just 0.81 G. Fig. [3] presents a statistical summary of our results, 
excluding measurements of non-HgMn stars. The upper panel 




1 2 

l<B z >/o(<B z >)l 




20 40 60 80 100 

a(<B z >), G 



Fig. 3. The upper panel: the distribution of the measured longitu- 
dinal field divided by its error. The lower panel: the distribution 
of longitudinal magnetic field error. 



shows the distribution of the ratio of (B z ) to its uncertainty. The 
lower panel illustrates the overall distribution of the longitudinal 
field uncertainties. The sensitivity of the majority of our mea- 
surements lies in the range of 1-20 G, which is much better than 
in previous studies of HgMn stars. 

For certain magnetic field configurations the first moment of 
the LSD Stokes V profile appearing in Eq. (|2]) may be zero, while 
the profile itself shows a magnetic signature. In order to exam- 
ine such situati ons we employed the False Alarm Probability 
(FAP) analysis dDonati et aU 119971) . It is a% 2 probability statis- 
tic, which estimates the probability that a given signal is pro- 
duced by random noise. We employed FAP analysis as a tool for 
the detection of a magnetic signal in the LSD Stokes V spec- 
tra, assuming no detection if FAP > 10~ 3 , marginal detection if 
10 3 <FAP< 10~ 5 and definite detection if FAP < 10~ 5 . We re- 
port no detection of any signatures indicative of a magnetic field 
in any of the 41 HgMn stars we observed. There is also no evi- 
dence of the coherent signal in any of the null spectra. For one 
star, HD 178065, the FAP analysis of the LSD Stokes V yielded 
a marginal detection. This result could be a statistical fluctua- 
tion since the measured FAP is close to the (arbitrary) boundary 
between no detection and marginal detection. No longitudinal 
magnetic field was detected for this star {{B z )= -2.26 + 2.08 G). 

A complete summary of our analysis of the circular polar- 
ization spectra of HgMn stars is presented in Table [2] In the 
columns from left to right we list the stellar HD number, des- 
ignation of a binary component (A - for primary, B - for sec- 
ondary and AB - for one measurement of the blended com- 
ponents of HD 27376). This is followed by columns contain- 
ing the Julian date of observation, the signal-to-noise ratio of 
the observed Stokes V spectrum and of the reconstructed LSD 
Stokes V profile, the mean longitudinal magnetic field derived 
from the Stokes V and null spectra, together with the respective 
uncertainties. In a few cases polarization measurements were 
based on 2 sub-exposures only, which did not allow us to com- 
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Table 2. Magnetic field analysis of HgMn stars. 



Target 


Comp. 


HJD-24 x 10 3 


S/N 


SVJV(LSD) 


<Sz)(V), G 


<S z )(null), G 


FAPxlO 


HD 1909 


A 


55202.527489 


200 


1494 


-4.02 ± 7.25 


9.64 ±7.11 


2.1790 


HD 1909 


B 


55202.527489 


200 


1494 


137.13 ±58.90 


-61.21 ±57.77 


0.4248 


HD 1909 


A 


55210.520641 


178 


1262 


4.60 ± 8.25 


-8.42 ±8.15 


1.0770 


HD 1909 


B 


55210.520641 


178 


1262 


47.95 ±61.14 


-99.79 ± 60.43 


9.6940 


HD 1909 


A 


55212.522266 


174 


1330 


16.81 ±8.87 


2.77 ± 8.76 


4.2570 


HD 1909 


B 


55212.522266 


174 


1330 


12.29 ± 56.73 


-48.85 ± 56.08 


5.6790 


HD 27376 


AB 


55201.768664 


507 


3355 


7.91 ±3.96 


6.82 ± 3.94 


4.5860 


HD 27376 


A 


55210.756062 


422 


2926 


-0.30 ±5.26 


10.76 ± 5.23 


4.0610 


HD 27376 


B 


55210.756062 


422 


2926 


-8.80 ±7.01 


-5.19 ±6.98 


2.2240 


HD 27376 


A 


55212.741449 


589 


4083 


-2.45 ± 3.89 


-2.22 ± 3.88 


9.3210 


HD 27376 


B 


55212.741449 


589 


4083 


4.78 ± 4.41 


3.55 ±4.40 


9.9540 


HD 27376 


A 


55213.730549 


411 


3071 


5.62 ± 4.91 


8.78 ± 4.90 


8.4040 


HD 27376 


B 


55213.730549 


411 


3071 


1.13 ±6.48 


-12.67 ±6.46 


3.9460 


HD 28217 




55213.538340 


283 


1472 


-12.30 ±46.87 


46.35 ± 46.60 


1.4140 


HD 29589 




55201.637379 


267 


1302 


-32.34 ±41.38 


-38.16 ± 41.16 


0.5699 


HD 33647 




55204.752546 


177 


1129 


29.51 ±38.47 


-34.16 ±37.99 


3.8950 


HD 33904 




55204.767578 


600 


4059 


-2.55 ± 2.59 


1.94 ± 2.57 


4.0470 


HD 34880 




55201.786748 


154 


933 


55.57 ±49.14 


79.04 ± 47.99 


9.2710 


HD 35548 




55202.740129 


172 


1588 


-0.74 ± 2.54 


0.29 ±2.51 


5.0180 


HD 36881 




55200.768514 


262 


2451 


-7.09 ± 5.49 


4.58 ± 5.47 


8.1070 


HD 42657 




55202.757839 


185 


1220 


-30.84 ± 79.97 


17.52 ±79.25 


7.4910 


HD 53244 




55204.779860 


477 


2692 


14.71 ± 13.06 


16.05 ± 12.98 


0.0655 


HD 53929 




55200.792070 


319 


1783 


7.14 ± 10.77 


-13.22 ± 10.59 


9.9250 


HD 63975 




55204.793412 


299 


1828 


8.28 ± 11.55 


-3.83 ± 11.34 


6.8330 


HD 65950 




55206.754842 


143 


1117 


-0.69 ± 16.87 


-24.08 ± 16.35 


7.2770 


HD 65950 




55211.765235 


152 


1039 


6.78 ± 18.33 


12.50 ± 17.77 


7.1000 


HD 65950 




55212.797482 


174 


1355 


-2.35 ± 13.79 


9.05 ± 13.51 


3.4760 


HD 68099 




55202.773866 


164 


1096 


-108.34 ±69.86 


40.66 ±68.10 


4.9430 


HD 68099 




55212.821947 


206 


1381 


96.70 ± 58.42 


31.10 ±57.31 


6.1680 


HD 70235 




55205.769179 


178 


1292 


22.95 ± 14.90 


2.99 ± 14.66 


4.9740 


HD 70235 




55212.837161 


204 


1581 


30.27 ±11.93 


-14.35 ± 11.72 


9.1010 


HD 71066 




55205.786739 


389 


2730 


-1.14 + 0.81 


0.24 ± 0.80 


0.6773 


HD 71833 




55205.812637 


193 


1253 


27.38 ± 26.46 


-14.94 ±25.91 


5.7850 


HD 71833 




55212.855644 


153 


1171 


-1.69 ±27.67 


20.61 ± 27.06 


8.1900 


HD 72208 




55202.789124 


140 


1289 


61.39 ±38.66 


-48.43 ± 38.20 


7.2000 


HD 72208 




55209.769108 


157 


1445 


-2.11 ±41.08 


-46.09 ± 40.41 


0.1644 


HD 75333 




55206.823570 


192 


2155 


4.38 ± 15.96 


-3.40 ± 15.78 


3.8920 


HD 78316 




55202.803837 


230 


1311 


0.03 ±4.12 


0.42 ± 4.04 


3.9410 


HD 78316 




55211.787325 


339 


1926 


-0.38 ± 2.66 


-3.40 ± 2.64 


5.3000 


HD 90264 


A 


55201.850566 


255 


1461 


3.43 ± 8.39 


-33.70 ± 8.23 


5.7250 


HD 90264 


B 


55201.850566 


255 


1461 


3.34 ± 12.42 


-2.50 ± 12.17 


0.6511 


HD 101189 




54982.974439 


51 


262 


-95.22 ± 65.63 


-53.75 ± 62.58 


0.4069 


HD 101189 




55201.862940 


179 


2797 


-6.41 ± 5.43 


4.52 ± 5.37 


1.4030 


HD 106625 




54983.067376 


149 


2493 


-15.99 ± 18.20 


13.06 ± 18.41 


4.1200 


HD 106625 




55205.872226 


881 


4348 


14.71 ± 11.11 


7.98 ± 11.13 


10.000 


HD 110073 




55205.856731 


417 


2813 


2.96 ± 4.25 


1.71 ±4.26 


9.9040 


HD 141556 


A 


54982.163224 


120 


1613 


-1.91 ±3.46 




7.4680 


HD 141556 


B 


54982.163224 


120 


1613 


8.38 ± 10.81 




1.4170 


HD 141556 


A 


55317.770736 


143 


1710 


1.54 ±4.48 


-12.63 ±4.51 


2.9260 


HD 141556 


B 


55317.770736 


143 


1710 


17.90 ± 11.08 


8.03 ± 11.16 


1.9210 


HD 149121 




55316.735917 


217 


1840 


-0.02 ± 4.45 


1.04 ±4.48 


6.4580 


HD 158704 




55316.845983 


172 


967 


-0.76 ± 3.56 


-1.91 ±3.50 


2.6810 


HD 165493 


A 


55316.925608 


144 


796 


4.77 ± 5.56 


-2.10 ±5.51 


8.1200 


HD 165493 


B 


55316.925608 


144 


796 


-31.05 ±27.67 


-80.86 ± 27.42 


0.9936 


HD 175640 




55317.865578 


170 


1327 


-0.59 ±2.16 


-3.05 ±2.16 


8.4280 


HD 178065 




55317.879714 


57 


431 


-1.92 ±8.81 




9.6770 


HD 178065 




55319.831679 


147 


1323 


-2.26 ± 2.08 


-4.12 ±2.08 


0.0085 


HD 179761 




55319.926735 


123 


1048 


23.92 ± 19.26 


-21.95 ± 19.07 


5.1960 


HD 186122 




55319.850529 


203 


1252 


0.66 ± 1.85 


0.45 ± 1.83 


7.3990 


HD 191110 


A 


55319.934190 


143 


1256 


1.67 ±7.17 


4.27 ±7.13 


9.8480 


HD 191110 


B 


55319.934190 


143 


1256 


6.26 ±9.15 


20.55 ±9.10 


7.2040 


HD 193452 




55316.939679 


90 


964 


-2.10 ±2.37 




2.2220 


HD 193452 




55319.868330 


203 


1868 


-1.04 ± 1.42 


-0.70 ± 1.41 


2.0460 


HD 194783 




54984.354843 


160 


1316 


10.08 ± 19.75 


-3.23 ± 20.05 


8.1140 


HD 194783 




55319.945379 


136 


1064 


44.01 ±21.60 


24.02 ±21.07 


4.0580 


HD 202671 




54984.382173 


175 


1220 


25.51 ± 16.01 


-25.08 ± 16.23 


2.3420 


HD 211838 




54984.396740 


150 


1432 


-194.23 ±77.70 


-2.01 ± 79.02 


7.1810 


HD 221507 




55210.548529 


349 


2498 


4.21 ± 7.04 


-9.51 ±6.99 


9.4710 
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pute the null spectrum. The uncertainties correspond to lcr con- 
fidence levels. The last column gives the FAP for the presence of 
a signal in the LSD Stokes V profile. All LSD Stokes /, V and 
null spectra are presented in Fig. [5] (Online material). 

4. Discussion 

As already mentioned, previous studies of magnetic fields in 
HgMn stars did not detect longitudinal fields stronger than 29- 
100 G, and in one particular case the limits were 6-19 G. These 
1 cr error bars are not low enough to exclude the presence of the 
longitudinal fields weaker than ~100 G, which could result from 
dynamically important local magnetic field structures covering 
only a fraction of the stellar surface. Thus, our understanding 
of the magnetic characteristics of these stars was fundamentally 
incomplete. In order to achieve a more detailed picture of the 
magnetic properties of HgMn stars, we studied a large sample of 
such stars with the most stable and modern spectropolarimeter 
available. 

Our analysis of the Stokes V spectra showed no evidence of 
any polarimetric signal, setting 3<x upper limits on the longitu- 
dinal field present of 3-30 G for the majority of the stars we 
observed. For one sharp-lined HgMn star we obtained an un- 
certainty of cr B , <1 G, which is better than previously achieved 
with a single observation of any early-type star. In addition to 
obtaining null results for longitudinal magnetic field measure- 
ments, we analysed the Stokes V LSD profiles themselves, find- 
ing no evidence for complex polarization signatures. Thus, we 
conclude that mercury-manganese stars do not possess strong, 
but topologically complex, magnetic fields. If we restrict our- 
selves to a dipolar magnetic field topologies, the dipole strength 
can be roughly estimated as B& ~ 3 (B z ). Based on our longitu- 
dinal field measurements, we conclude that the HgMn stars do 
not possess dipolar fields stronger than 3-30 G. 

It is important to evaluate our results in the context of 
the possible role of magnetic fields in the recently discovered 
spot formation on HgMn stars. Although there are no detailed 
models describing the process of the interaction bet ween mag- 
netic fields and the plasm a, it is understood (e.g., IWade et aD, 
20061 lAuriere et al., 2007) that this interaction becomes impor- 
tant when the strength of magnetic field exceeds the equiparti- 
tion limit, defined by I Wade etaf] d2006l) as B eq = ^\2nP gas . At 
this field strength the energy of the magnetic field equals that of 
the gas. 

We estimated the value of the equipartition magnetic field 
strength for atmosphere models with effective temperatures of 
T s g= 10000 K and r e ff= 15000 K, representing the range covered 
by our sample of HgMn stars. Using thermodynamic parame- 
ters tabulated in these models and assuming the ideal gas law, 
we calculated the equipartition field as a function of the optical 
depth in the continuum. The result is presented in Fig. [4] From 
this figure we can see that the value of B eq at log T5000 = -0.5 
is ~200 G for the stars of both low and high effective tempera- 
tures. The value of B e<i then decreases to 100 G at log T5000 = -2 
for T eS =10000 K and to 60 G for T eff =15000 K. The majority 
of spectral lines used in our measurements form over a range in 
optical depth from log T5000 = -0.5 to -2. In this interval the 
values of B eq are much higher than our upper limits for the or- 
ganized magnetic field on the stellar surface. Thus, it appears 
there are no magnetic fields on the surfaces of HgMn stars that 
could lead to the formation of chemical spots. These conclusions 
apply to the fields stronger than ~ 10 G probed by our study. On 
the other hand, we cannot exclude the presence of complex mag- 
netic fields with a field strength of ~ 1 G in surface spots, similar 




I I I I I I I 1 



-6 -5 -4 -3 -2 -1 

log ^5000 

Fig. 4. The equipartition magnetic field as a function of the op- 
tical depth. The solid line corresponds to the cooler model and 
the dashed one to the hotter. 



to the tangled magnetic topology suggested for Vega (Peti t et al.L 
1201 01) . 

These results lead to the conclusion that there exists a real 
dichotomy between HgMn and magnetic Ap stars. The latter 
show chemical inhomogeneities in the presence of dipolar-like 
magnetic fields stronger than 300 G. This limit is determined by 
the field strength that can withstand shearing by di fferential rota- 
tion a nd survive at the surfaces of early-type stars dAuriere et al.L 
120071) . In contrast to magnetic Ap stars and despite a similarity 
in fundamental parameters, HgMn stars clearly do not have mag- 
netic fields reaching this limit. Consequently, the surface chem- 
ical structures in mercury-manganese stars must form under a 
mechanism fundamentally different from magnetically-driven 
spot formation in hot Ap or cool active stars. 

A detailed Doppler imaging study of the surface struc- 
tures on the brightest HgMn star, a And, carried out by 
Koc hukhov et a l. (2007) using observations collected over a pe- 
riod of seven years led to the remarkable discovery of Hg spot 
evolution. This is very different from the spotted early-type mag- 
netic stars, in which the morphology of spots remains steady 
over te ns of years. The fa ct that no magnetic field was detected in 
a And (IWade et all 12006) makes the explanati on of the spot for- 
mation and evolution a very challenging task. iKochukhov et all 
(I2007h suggested that mercury atoms are concentrated in a thin 
layer where the radiation pressure and gravitational force com- 
pensate for each other. The observed changes in the spot mor- 
phology are then ascribed to hydrodynamical instabilities, upset- 
ting the fine force balance in the Hg cloud. The physics of this 
intricate process and its relation to the binarity, stellar rotation 
and dynamic atomic diffusion effects is yet to be explored. 

Despite being frequent targets of detailed abundance analy- 
ses during more than half a century, HgMn stars were recognized 
as spotted variables only recently. To date, 7 spotted HgMn stars 
have been found. It is uncertain how many more such stars will 
be discovered. Prior to the formulation of exotic theories of non- 
magnetic chemical spot formation, it is important to determine 
a complete picture of the occurrence rate of chemical inhomo- 
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geneities on the surface of HgMn stars. If spots are found in 
only a small fraction of these objects, it may be that the ob- 
served structure formation is associated with a certain stage of 
stellar evolution or a certain configuration of the binary system. 
On the other hand, discovery of spots in the majority of HgMn 
stars would indicate the presence of a universal hydrodynamical 
process, not anticipated in current stellar models. 

We believe that in-depth studies of selected spotted HgMn 
stars as well as a broad search for spectral line variability in 
all accessible objects of this type will provide sufficient obser- 
vational constraints for the development of a theory of non- 
magnetic spot formation. Analysis of known spotted stars can 
yield information on the time-scales of spot evolution and, pos- 
sibly, on the heights of these inhomogeneities relative to the typi- 
cal line formation regions. It is also important to continue search- 
ing for weak magnetic fields in spotted HgMn s tars as previ- 
ously done bv lWade et al.1 (120061) for a And and bv lFolsom et alJ 
(2010) for AR Aur. On the other hand, a line profile variability 
survey of a large sample of HgMn stars will help us to uncover 
dependencies on stellar parameters, rotation rate or binarity, pos- 
sibly revealing a correlation of one these factors with the pres- 
ence of spots. In subsequent studies of HgMn stars we plan to 
carry out both detailed studies of individual objects and to per- 
form variability surveys to find new spectrum variable stars. 
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Fig. 5. From bottom to top each panel shows the LSD profiles of Stokes I, V and null spectrum (if available) for all HgMn stars in 
our survey. 
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Fig. 5. Continued. 
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Fig. 5. Continued. 
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Fig. 5. Continued. 



